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Abstract In this work we review the efforts that have been done to study the exci- 
tation of the standing fast kink body mode in coronal loops. We mainly focus on the 
time-dependent problem, which is appropriate to describe flare or CME induced kink 
oscillations. The analytical and numerical studies in slab and cylindrical loop geome- 
tries are reviewed. We discuss the results from very simple one-dimensional models to 
more realistic (but still simple) loop configurations. We emphasise how the results of the 
initial value problem complement the eigenmode calculations. The possible damping 
mechanisms of the kink oscillations are also discussed. 
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1 Introduction 

Transversal coronal loop oscillations are regularly observed by the EUV telescope on 
board TRACE. Nowadays, there is abundant information about the period, d amp- 
ing time and amplit ude of the oscillations (see for example lAschwanden et al. 2 003 : 
ISchrii ver et al ] |2002l . and see also the review of M. Aschwanden in this volume) . In most 
of the cases these oscillations are triggered by a flaxe or a CME, and they have been in- 
terpreted as fast standing magnetohydrodynamic (MHD) kink oscillations of cylindrical 
tubes (see INakariakov and Ofmanll200]] : iRuderman and Robert^l2002l : [Coossens et al.l 
I2OO2I '). 

In general, most of the theoretical studies on loop oscillations assume that the 
system is in the stationary state and that the loop dynamics is given by the normal 
modes. The calculation of the eigenmodes is possible in rather simple configurations but 
their properties are usually unknown in complex structures. Although normal modes 
should be seen as the building blocks to interpret coronal loop oscillations they do 
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not represent the whole picture, but their study provides a basis for understanding 
the dynamics of the system. To have a more accurate description, the time-dependent 
problem needs to be analysed. Ideally, we would like to know which normal modes are 
excited given an initial perturbation. 

Although there is an extensive literature about MHD oscillations in coronal loops 
only those theoretical works where the time-dependent analysis has been considered 
are reviewed here. Usually the time-dependent problems are classified in driven or in 
initial- value problems. Driven problems are useful to model, among others, the effect of 
the photospheric motions in coronal loops. In this case a continuous motion is applied 
at the base of the loop exciting several kinds of modes that propagate along the tube. 
On the other hand, initial-value problems are more appropriate to describe a sudden 
release of energy, as in flares. Such energy release generates a disturbance, a blast wave 
or an EIT wave, that propagates in the corona and produces kink oscillations of nearby 
loops. In this work we will concentrate on initial-value problems. Moreover, since most 
of the observations of kink oscillations are interpreted as the fundamental standing kink 
mode, here we will focus only on this mode, characterised by a single frequency and a 
dominant wavenumber along the loop. For propagating kink modes, characterised by 
a wide superposition of different longit udinal wav eimmbcrs and fr equencies, the reader 
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is referred to the works of.Murawski and Robertsi (,1993a. h , 19941') • Hua.ng et al.l (Il999l l 



l|2004h : iTerra-Homem et"al] l|2003l ) 



and to 
cylindrical loops. 

In this work we flrst describe the results of the time-dependent problem starting 
with simple loop models and then we increase the complexity of the equilibrium. We 
also centre on those studies where a direct comparison with eigenmodes can be done 
since this helps to understand and interpret the results of the time-dependent problems. 



2 Basic Loop Models 

Two basic equilibrium models have been studied in detail in the literature. They are 
the simplest representation of coronal loops and are important to understand the basic 
properties of MHD waves. These models are the slab and the cylindrical flux tube. 



2.1 Slab 



The simplest loop model is a Cartesian slab configuration. In this model the loop 
is represented by a plasma column of enhanced density, p;, respect to the coronal 
environment, with density pc. The slab has a half- width of a and the magnetic field is 
vertical (B = Bgez). Since the inhomogeneity is in one direction, hereafter assumed in 
a;, it is possible to make Fourier analysis in t he ot her i gnorable directions {y a nd z). 
The reader is referred to lEdwiri and Robert^ (|l982l ) and lTerradas et al.l (|2005al ) (with 
special attention to the leaky modes) for the details about the dispersion diagram in 
this model. When the wave propagation is in the plane of the slab then the frequency 
of the fundamental fast kink body mode, in the limit of thin tubes (when k^a ^ 1), 
tends to the external Alfven frequency. 



(1) 
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where uac ~ Bq/^hPo is the external Alfven speed. However, when the perpendicular 
pro pagation is allowed (ky ^ 0) the body kink m ode is transformed in a surface wave 
(see IZhelvazkov et all 1 19961 : lArregui et al] |2007| ) and its frequency tends to the kink 
frequency (when ky ^ kz) 



.,^k..b^^(^. (2) 

V Pi + Po 

In general this model is appropriate if the structure is invariant in the perpendicular 
plane (i.e. in the j/— direction). 



2.2 Cylinder 

A cylindrical f lux tube is the natura l extens i on of the p revious model. The reader is 
referred to Edwin and RobertsI (|l983l ): ICallvl (|l986l . liooi l for the details of the disper- 



sion diagram and the classification of the different modes (trapped and leaky). Here 
we will restrict to body waves and the m = 1 mode (the oscillation that displaces 
transversally the whole tube). In the limit of thin tube the frequency of oscillation is 
the kink frequency, given by equation ([2|. 

The main advantage of the cylindrical model is that it represents a three-dimensional 
structure. It is interesting to note that the kink body modes are very similar to the 
kink modes of the slab model when ky is large. The frequency of oscillation is the same, 
the kink frequency, and the eigenfunc tion is quite localised ne ar the loop boundary. In 
fact in some papers (see for example iHollweg and Yang||l988l ). the results of the slab 
models have been extended to the cylindrical geometry by using the equivalent of the 
azimuthal wavenumber in Cartesian coordinates, i.e. assuming that ky = m/R. 

On the other hand, one of the assumptions that is usually relaxed in the cylindrical 
and in the slab model is the Alfven speed profile across the boundary. In the simple 
models it is assumed that the density (and the Alfven speed) has a discontinuous 
jump at the interface between the loop and the coronal environment. However, if the 
Alfven speed has an inhomogeneous layer, with a typical width of I, then the process 
of resonant absorption takes place (for m > 1 in the cylindrical model or for ky ^ in 
the slab model). 



2.3 Curved configurations 

One of the configurations that takes into account the curvature, neglected in the previ- 
ous models is the circular loop model. In this force-free model the shape of the magnetic 
field lines is given by a purely poloidal magnetic field, 

Se=Bo^. (3) 
r 

The density profile is usually chosen (in the zero-/? approximation) to have the same 
dependence as the magnetic field and with a density enhancement, representing the 
loop, at a given radius. The advantages of this model is that it is possible to make 
Fourier analysis along the direction and the loop cross section is constant along the 
tube. One of the disadvantages is that it has a singularity at r = and this might be 
difficult to handle from the numerical point of view. 
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Another popular model is the potential arcade. In this model the current-free mag- 
netic field is given by 

B:, = Bosin(x/ylB) e-^/^«, (4) 
Bz = -Bo cos (x/Ab) e""^/-^^ , (5) 

where Ag = L/n and L is the arcade width. The background density is usually taken 
to be of the form 

(6) 

The parameter A allows to have a decreasing or increasing Alfven speed with height (v^^ 
is independent of the a;— coordinate) and the loop is modelled by a density enhancement 
along certain field lines. The configuration is well behaved at the origin but the loop 
cross section is not constant along the tube and no Fourier analysis is allowed along 
the field lines. 

These two equilibrium models are two-dimensional but they can still be used to 
study curved loops in three dimensions assuming that the field is invariant in the 
perpendicular direction to the plane. Other sophisticated equilibrium models involve 
three-dimensional dipoles or even photospheric magnetic field extrapolations, usually 
force free, together with the inclusion of a certain density profile along the field lines. 
However, these models are still in a very primitive stage. 



-z A 



3 Waves in Slabs 

3.1 Basic properties of fast MHD waves in a uniform line-tied medium 

We start with the most simple configuration, a uniform plasma permeated by a vertical 
and uniform magnetic field. In the P zero limit the ideal linearised MHD equations 
reduce to the well-known wave equation for the fast modes, 

where ^ is the horizontal displacement and vp^ = BqI ^fjip is the Alfven speed. Since 
/3 = slow modes and Alfven waves (we assume that ky — 0) are absent. We Fourier 
analyse in the z-direction, i.e. we assume that the perturbations are of the form e~''^^^. 
Hence, equation ([71 reduces to the Klein-Gordon equation, 

where uc — vp^kz is the cut-off frequency. Thus, after introducing the longitudinal 
wavenumber kz the two-dimensional wave equation reduces to a one-dimensional equa- 
tion, and by selecting the appropriate value of this parameter the effect of line-tying 
can be easily incorporated. For the fundamental oscillation with a node at each foot- 
point we have that kz = it/ L, L being the length of the system (or loop if we include a 
density enhancement) in the vertical direction. In this model the time a wave needs to 
travel along the entire loop to feel the effect of the photosphere is neglected, and the 
loop instantly feels the effect of line-tying. 



(7) 



(8) 



5 



Normal mode solutions of equation ([8]), i.e. with a dependence of the form (.^i^^^^^) ^ 
yield the dispersion relation lo = ^J^v^p^ + which indicates that wave propagation is 
dispersive. The propagation speed for waves with wavenumber k is the group velocity, 



which ranges between zero, for waves with infinite wavelength, and wa for waves with 
infinitely small wavelength. Therefore, modes with short wavelengths propagate faster 
than modes with larger wavelengths. In addition, any mode with frequency below the 
cut-off frequency is evanescent and unable to propagate in the medium. The dispersive 
behaviour of waves is simply due to the fact that we assume a fixed wavelength along 
the loop (a single kz) to model the effect of fixed footpoints. 
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Fig. 1 Displacement as a function of time at position x = 2L for different wave source widths: 
the dashed line corresponds to a (5-function and the continuous line corresponds to a Gaussian 
profile with w = O.IL. The numerical solution with a Gaussian profile with w = O.IL for a thin, 
low density slab model (a = O.OIL and Pi/pe = 2) is displayed with circles. For comparison 
purposes all the profiles have been normalised to the same amplitude. It is clear that different 
kinds of impulsive excitation and the presence or not of the loop result in essentially the same 
temporal behaviour for this low density contrast, thin loop. 




Once we know the properties of fast MHD waves in the model we can study the 
time-dependent problem by launching a perturbation in the system. An analytical 
solution of equation (|8]) can be derived for the following initial conditions 

e(a;,0)=0, ^{x,0)^vo5{x), (10) 
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vq being the velocity of the initial perturbation. The displacement as a function of time 
in terms of the Bessel function Jg is given by 

a^,t) = £^Jo {u^cVt^ - i^/vAr) m^xM, (11) 

where TL is the Heaviside function. This solution (see for example also lRae and Roberts! 
Il982l ). is characterised by a wave front travelling at the Alfven speed in the direction 
normal to the magnetic field lines and by a wake behind this pulse oscillating at the 
cut-off frequency Uc (see dashed line in Fig.[T}. For t ^ x/i^a the asymptotic expression 
for the displacement is 

t) = ^ "° 2 , cos(c^ct - n/A), (12) 



1/9 

which indicates that the displacement amplitude decays as t ' and the frequency 
of oscillation is simply ujc- Note that the atten uation has noth i ng to do with any 
dissipation mechanism. These re sults, derived by iTerradas et all 1 2005bl ) are slightly 



different from those obtained bv lUralovl ( 20031 ) . The solution presented here 



IS more 



appropriate to describe the physical problem of a fiare-induced perturbation and the 
propagation of the resulting disturbance in the x-direction. 

In order to have simple analytical expressions we have considered that the initial 
perturbation is a (5-function, but it is possible to consider a source of finite width (for 
simplicity we assume a perturbation with a Gaussian profile of width w). In Figure [T] 
we have represented the solution for w — O.IL (see continuous line). For a typical 
flare event the spatial scale of the compact flare source is much smaller than the mean 
length of the oscillating loops, hence values of w below O.IL are close to realistic source 
widths. For large t the solution tends to the results of the (5-function perturbation and 
the differences are mainly in the wave front but not in the oscillating wake. Hence, 
an excitation with a Gaussian shape does not change much (for a reasonable range 
of disturbance widths w) the displacement profile in comparison with the 5-function 
perturbation. 



3.2 Effect of structuring: Leaky and trapped waves 

Now we include a density enhancement (of half-width o) representing the loop. In this 
slab model, the Alfven speed inside the loop is smaller than the Alfven speed in the 
corona (vai < fAo)- Since the Alfven speed changes in the x-direction the analytical 
solution of the initial value problem , involving temporal Laplace transforms (see for 
example I Andries and Goossenj [20071 ) . is much more complex. For this reason we hav e 
solved equation ^ numerically using standard techniques (see ITerradas et "al]|2005al ). 
To excite the kink mode we assume the following perturbation 



Vx{x,t = 0) = Vxo exp 



(13) 



where we choose w = a (now instead of the displacement, ^, we use the velocity, Vx 
(vx = ^). The numerical solution (see Fig. [2]) shows that the induced disturbances 
propagate away from the slab and contain different wavelengths due to the dispersive 
nature of fast waves described before. However, the shape of the velocity inside the slab 
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and its near surroundings has a different behaviour due to the presence of the cavity 
or wave guide. It resembles the form of the fundamental kink mode eigenfunction: vx 
has an extremum at a: = and decreases exponentially outside the loop. This indicates 
that the slab is oscillating with the normal mode. 
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Fig. 2 Plot of for different times. The initial perturbation is given by equation JTSj with 
•w = a and a = 0.05L. Note the different spatial scale for t = 0. The grey area represents the 
density enhancement (pi/po = 3). 



In Figure [3] we have plotted Vx as a function of time at the slab centre (a; = 0) . 
The frequency of oscillation for large times is in excellent agreement with that of the 
fundamental kink mode calculated by solving the dispersion relation of the slab (it 
is basically given by eq. [1]). We also see that there is a short time interval before 
the loop reaches the stationary regime. In Figure |3] we have plotted a detail of this 
transient, which is initially dominated by the propagation of the initial disturbance 
inside the slab (0 < t/rA < iajvp/) and is followed by the excitation of a leaky mode 
(for 2a/vjs^i < t/rj/j because of the reflection of the disturbance at the slab boundaries. 
During this impulsive leaky phase the amplitude shows short period oscillations. Part of 
the initially deposited energy in the slab is radiated away through the fundamental kink 
leaky mode before the loop oscillates with the corresponding trapped mode. The period 
of these oscillations, which are now superimposed to the much large r period of the 
trapped kink mode, agrees with the period of the first leaky mode (see iTerradas et ahl 
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Fig. 3 Plot of Vx at the slab centre (x = 0) for the numerical simulation of Figure [2] After a 
short transient phase the loop reaches a stationary state, oscillating at the frequency predicted 
by the normal mode analysis. In the inner plot the short impulsive leaky phase is displayed, 
showing the signal for t/rA between and 1. 



given by the following equation (in the limit a/L <Si 1) 



2a 1 

VAi n 

2a 



In 



1 + ^Ai/wAo 



1 - WAi/wAo 



(14) 
(15) 



Note that the damping time {t^ = lo^^) is independent of n, in the limit a/L <^ 1 all 
leaky harmonics (i.e. modes with different n) have the same damping time. Moreover, 
when t^Ai ^ I'Ao the damping time reduces to 



Td-a-j-. (16) 
"Ai 

The time-signatures in the slab model are now patent, after an initial disturbance 
there is a short leaky transient followed by the excitation of the trapped mode, which 
oscillates with constant amplitude with time (in ideal MHD, i.e. without any dissipation 
mechanism) . 

It is worth noticing (see Fig. [3]) that the amplitude of oscillation of the normal 
mode is constant but smaller than the initial amplitude (vxo)- K the perturbation is 
located in the external medium much less energy is trapped in the cavity. In this case 
the motion of the slab is dominated by the wave passing through it (described in the 
previous section) rather than by the eigenmode. In Figure[T](see circles) we have plotted 
the evolution of the slab for a perturbation located in the corona. If we compare this 
result with the profile when there is no density enhancement (see continuous line) we 
find that the curves are almost the same. Thus, for an external excitation the slab is 
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Fig. 4 Plot of Vx for different times. The initial perturbation is given by equation JTSj with 
w = 2a and a = 0.02L. In this simulation kya = 0.5. Note the different spatial scale for i = 0. 
The grey area represents the density enhancement (pi/pe = 10). 



unable to trap a significant amount of energy. This issue is discussed later in detail for 
the cylindrical loop model (Section 14. 2|l . 

The effect of oblique pro pagation (ky ^ 0) in t he slab, not included in the previous 
models, was investigated by lArregui et al. I l|2007l ). The properties of the eigenmodes 
in this system are interesting since for large ky the body kink mode is converted into 
a surface kink mode. This is clear in the time-dependent solution plotted in Figure |4] 
for kya = 0.5. The system is oscillating with a trapped mode but the surface nature 
of the mode is revealed in Figure |4}i, where we can appreciate that the amplitude 
of the oscillation is highly localised around the loop boundaries. Moreover, the fre- 
quency of the oscillations tends to the kink frequency (given by eq. [2]). The effect of 
oblique propagation has been extensively studied in the past specially in the context 
of resonant absorption a t magnetic interfaces (see for example iLee and Robert^ll986l : 
iHoUweg and Yang||l988l ). 



3.3 Multi-structures 

Many loops in the solar corona are not isolated but forming groups or bundles of loops. 
For this reason we now consider an idealised system of two parallel loops modelled as 
dense plasma slabs of half-width a and length L. The distance bet ween the centres of 
the slabs is d. Such a configuration was studied bv lLuna et al. The normal mode 
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analysis indicates that there are two kinds of normal modes (see also lPfaz et al.l[2005l ): 
solutions symmetric with respect to a; = (both slabs move in phase) and solutions 
antisymmetric with respect to 2; = (the slabs move in anti-phase). The symmetric 
mode respect to the centre of the structure is the only trapped mode for all distances 
between the slabs while the antisymmetric mode is leaky for small slab separations. 
Nevertheless, there is a wide range of slab separations for which the fundamental sym- 
metric and antisymmetric trapped modes are allowed and have very close frequencies. 
These modes are excited according to the parity of the initial perturbation. 
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Fig. 5 Time-evolution of Vx for d = 0.5L and a non-symmetric initial excitation. The grey 
areas represent the density enhancements (pi/pc = 3). 



In Figure [5] the results of an initial excitation for a loop separation d = 0.5L are 
represented. The initial disturbance is half the sum of the symmetric and antisymmetric 
initial conditions. During the initial stages of the temporal evolution (Figs. [5}d and c) Vx 
has no definite symmetry with respect to x = because the solution is the sum of the 
symmetric and antisymmetric modes. These modes are the fundamental symmetric and 
the fundamental antisymmetric, which are, for this particular slab separation trapped 
and leaky, respectively. As a consequence, after some time (Fig. [5}1) the antisymmetric 
mode amplitude is negligible in the vicinity of the slabs and the system oscillates 
in a symmetric manner. The oscillation is collective and different from the motion 
associated to the individual eigenmodes of the slabs. Rather different results are found 
when the system is excited with the same initial condition but for a distance between 
the slab centres of d = 2L. This choice of the slab separation results in the fundamental 
antisymmetric mode becoming trapped. The evolution of the system is plotted for 
different times in Figure [6] and, although after some time the two slabs seem to move 
in phase (see Fig. (6)3) , in a later stage the right slab has transferred all its energy to the 
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Fig. 6 Time evolution of for d = 2L and a non-symmetric initial disturbance. Note the 
interchange of energy between the two slabs in the last two frames. 



left slab and is motionless (see Fig.|6};). At an even later time (see Fig.[6}i) the picture is 
just the opposite, with the left slab fixed and the right slab in motion, this is the result of 
the continuous exchange of energy. In Figure [7^ we have represented Vx at the centre 
of both slabs. Contrary to the behaviour in the stationary regime for symmetric or 
antisymmetric initial perturbations, the oscillations do not attain a constant amplitude, 
but they instead display a sinusoidal modulation. This is a beating phenomenon due 
to the simultaneous excitation of the symmetric and antisymmetric modes with similar 
frequencies. These frequencies are recovered from the power spectrum of the velocity at 
the centre of right slab (see Fig.[7}3), which shows two power peaks with periods almost 
identical to those of the fundamental antisymmetric eigenmode and the fundamental 
symmetric eigenmode. 



3.4 Curved configurations 

One of the main properties of the straight slab models described before is that the 
system always has trapped modes. However, when the configuration is more complex 
then the behaviour of the solution can be leaky. For example, in two-dimensional curved 
configurations the Alfven frequency changes (in general) with height, and the phase 
speed of the body mode might be above the local Alfven speed at a certain position 
in the arcade. The character of the wave is propagating and thus the mode is leaky, 
being its amplitude attenuated with time due to the continuous radiation of energy. 
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Fig. 7 (a) Vx measured at the centre of the slabs for the simulation shown in Figure [6] The 
solid and dashed lines correspond to the right and left slabs, respectively, (b) Power spectrum 
of the previous signal (solid line). The two vertical dashed lines indicate the periods of the 
fundamental antisymmetric mode and the fundamental symmetric mode calculated from the 
eigenvalue problem. 



3.4-1 Circular slab 



The first time-dependent analysis of the ex citation of (vertica l ) fast waves in a curved 
slab (described in Section [2.311 was done bv lBradv and Arbed l(2005l '). They considered 
a driven problem at one footpoint (and line- tying conditions at the other), and several 
fast body modes were excited. In Figure [8^ there is an example of the excitation of 
a particular longitudinal harmonic. We can appreciate that the energy is not only lo- 
calised around the slab (with footpoints around x = ±1), there are many fronts outside 
the slab that are a consequence of the leaky character of the mode. The a ttenuation of 
the oscillations with time due to the wave leakage is evident in Figure [8b. iBradv et al.l 
12006) extended the previous model and considered in more detail the effect of wave 
tunnelling on the ver tical fast modes. Later, t he eige nmodes of the circ u lar slab were an- 
alytical ly derived bv lVerwichte et aL I (|2006al lbl) and lOiaz et al. I l|2006l ).l Verwichte et al.l 
( 2006bh showed that the calculated dampin g rate of fast ma g netoa coustic waves is 
consistent with the numerical simulations of iBradv and ArbeJ ( 20051 ) . so a clear link 
between the eigenmodes and the time-dependent solution could be established. 




Fig. 8 (a) Velocity at the end of the simulation after the driver has stopped, (b) Decay of the 
velocity inside the slab as a function of time due to the emission of leaky waves. The system 
shows three distinct phases, the stable regions due to the driver, a rapid decay phase when the 
driver is switched of and the stable d e cay o f the normal model of the tube. These calculations 
were performed by IBradv and Arberl | |2005[) . 
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3.4-2 Slab embedded in a corona 



The main properties of propagating fast MHD waves in coronal arcades witho u t a den - 
sity enhancement have been described bv lOhver et"all (| 19981 ) : iTerradas et all (l2008bh . 
Recently, a battery of papers has been devoted to the analysis of the p roperties of these 
waves when a slab is embedded in the arcade. iMurawski et al.1 hoO^ ) showed that an 
impulsively disturbance near the apex of the loop is able to excite standing kink waves. 
Since the perturbation was very localised sev eral standing kink modes were excited at 
the same time. Similar results were found bv lSelwa et "aiT(|2005l ). who was able to ex- 
cite the fundamental standing kink mode with a pertur bation below the loop apex. 
A more detailed analysis was performed by ISelwa et al.l (|200dl who studied how the 
excitati on of the kin k mod es depends on the position of the pulse and its width. More 
recently ISelwa et al. I (|2003) have extended the previous works and have shown that in 
the arcade leaky waves are also excited. In Figure |9^ we can appreciate several wave 
fronts propagating upwards which are related to the emission of leaky waves (see also 
Fig.EK). It should be noted that these leaky waves are slightly different from those dis- 
cussed in Section [3. 2 1 However, the short period leaky transients should be also present 
in the curved case, but apparently there is no clear indication of such oscillations in 
the previous numerical studies. 

Interestingly, in all the simulations a strong damping of the vertical kink oscillations 
is always reported (see Fig. [IJa). Although in the arcade model there is leakage due 
to the process described in the previous section, it is not clear that it can account 
for the strong damping found in the simulations. Unfortunately, up to kn ow there are 
no av ailable eigenmode calculations of a slab in the potential arcade f see ISmith et al.l 
Il99i for some preliminary attempts) which would help to clarify the origin of the 
damping. It is also interesting to mention that in many simulations the loop does not 
return to its original position, as can be appreciated in Figure [9j3. This issue needs to 
be investigated in detail, since in the linear regime it is expected the loop to oscillate 
around the initial equilibrium. However, if the amplitude of the initial oscillation is 
large enough, the system can relax to a new magnetohydrostatic equilibrium. 
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Fig. 9 (a) Perturbed energy density. The position of the curved loop is mar ked with the 
black lines, (b) Density evolution at the loop apex. These results were derived bv lSelwa et al.l 



On the other han d, the influence o f the d ense photosphere on vertical oscillations 
has been analysed bv ICruszecki et aL I l|2008l ). By including a dense photosphere-like 
layer, instead of line-tying conditions, these authors claim that in this model there 
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is a more efficient excitation and attenuation of tfie vertical kin k mode due to the 
leakage through the photosphere. IGruszecki and Murawskil ( 20081 ) have analysed the 



effect of gravity on the oscillations and they have shown that the reported damping rate 
decreases as the coronal scale-height increases. However, this damping rate is still small 
compared with the observatio ns. On the other h and, t he effect of multi-structure in the 
arcade has been modelled by IGruszecki et al.l ( 200g) by including a double stranded 
loop and up to five strands. 

It is interesting to note that leakage in these two-dimensional structures is probably 
overestimated in comparison with three-dimensional loop models. In three-dimensions 
the eigenfunction of the kink mode is much more localised near the tube and the effect 
of inhomogeneities i n the external medium is presumably much less important (see 
iTerradas etlolbood ). 



4 Waves in Cylinders 

We turn now our attention to the analysis of kink body waves in the cylindrical loop 
model. 



4.1 Leaky and trapped waves 

The linearised MHD equations are solved in cylindrical coordinates, a nd since we are 
interes ted in the kink mode, we impose m = 1 in the equations (see ITerradas et al.l 
l2007bh . As in the slab model we perturb the system with an external localised distur- 
bance with the following form 



Vr{r, t = 0) = Vj.Q exp 



(17) 



where rg is the position of the centre of the disturbance and ui is a parameter related to 
its width. The results are qualitatively similar to those found in the slab configuration. 
The initial perturbation induces disturbances that propagate towards the loop and 
also wavefronts travelling in the opposite direction. Part of the energy of the initial 
perturbation is trapped while another part will be radiated through the excitation of 
the leaky modes. The radial velocity component at the centre of the loop is plotted in 
Figure [TOk . The signal initially shows a short transient phase (0 < t/r^ < 20) which 
is followed by a long-period oscillation. This oscillation is due to the excitation of the 
trapped kink mode. Its frequency is simply the kink frequency of the tube (see dashed 
line in Fig. llOb l. The amplitude of the oscillation of the trapped mode remains constant 
with time. On the other hand, the initial transitory phase is associated to the leaky 
modes, represented in more detail in Figure [TTk . The period of the transient phase is in 
agreement with the period of the (fundamen tal) kink leaky trig mode (see dashed line 
in Fig.int)- According to lCallvl (|l986l . booaj ). the trig mode has the following period 

(18) 



(rn + 2n+ |) i;a 



n being any integer consistent with P > 0. Their damping time is 



Td^pl^. (19) 
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The previous expressions are valid for VAi ^ ^Ac (and for thin and long loops) and 
are the modifications due to the cylindrical geometry of those derived for the slab (see 
eqns. [14] and [H]). Again for these leaky modes the period and damping time are 
independent of the loop length. Additionally, in the dispersion diagram there are other 
families of trig leaky modes, for example, the "Type II trig modes" and the "orphan 
mode" rsee lCallvl[l98i .[ 20031 ). However, these modes seem to be very difficult to excite 
by an initial perturbation. Something similar happens with the Principal Fundamental 
Leaky Mode, although it is a solution of the dispersion diagram (with a frequency very 
similar to that of the trapped kink mode for kzR <C 1) it is very difficult to excite by an 
initial disturbance since there is no clear evidence of such modes in the time-dependent 
solutions (see Tcrradas ot al. 2 007bl, for a detailed investigat ion about the excitation 
of this peculiar mode). In fact, iRuderman and Roberts! 1 20061 ) claim that this mode is 
spurious. 




150 200 



Fig. 10 a) Radial velocity as a function of time at r = 0. b) Wavelet transform of the signal. 
After a short transient phase, the loop oscillates with the frequency of the kink mode. The 
loop length is L = 30 /? and the density contrast is Pi/po = 3. The initial perturbation if given 
by equation I I17I I with ro = 5 R, w = R. The dashed line represents the frequency of the kink 
trapped mode calculated from the dispersion relation. 



4.2 Energy deposition by initial disturbances 

Once we known the time signatures of the an external disturbance on the loop, we 
turn to the theoretical study of how the amplitude of oscillation depends on the initial 
disturbance. This problem has been partially addressed in the slab section, where it was 
suggested that an external disturbance might deposit a very small amount of energy 
i n the loop. To analyse this problem we refer to the work of IRuderman and Roberta 
(|2006l ) who solved the initial value problem using the Laplace transform to determine 
the motion of the tube. These authors have found that the asymptotic behaviour of the 
loop for t ^ oo is given by the fundamental normal mode of the structure. Since we 
are mostly interested in the energy that is trapped in the loop and not in the transients 
before the loop settles in the normal mode (already discussed), instead of solving the 
full time-dependent problem numerically and c alculating the deposited energy in the 



loop we use an approach based on the results of !Ruderman and Roberta ( 2006! ) . These 
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Fig. 11 a) Detail of the radial velocity plotted in Figure [TOl in the range < t/rp^ < 20. 
b) Wavelet transform of the signal. The leaky mode is clearly identified and there is some 
indication of the trapped mode in the power spectrum. 



authors showed that the amplitude of oscillation of the loop is basically proportional 
to the con volution between the i nitial disturbance and the eigenmode of the loop (see 
eqn. (7) in lXerradas et al"]|2007ah . Using this method we can easily calculated how this 
a mplitude (i.e. th e energy of the trapped mode) depends on the different parameters. 
In lTerradas et al] ( 2007ai ) several types of perturbations (always in the total pressure) 



were studied. Since we are interested on the kink modes the simplest perturbation is 
an m = 1 disturbance localised in the external medium. 

Po^PNe-(^) , (20) 



where Pj\[ is a normalisation constant, rg is the location of the Gaussian and a is 
the width at half height (do not confuse with the loop radius R). This allows us to 
study the behaviour of the loop for different locations of the disturbance and different 
widths of the initial pulse. The profile of the initial perturbation for m — 1 using the 
radial dependence given by equation (|2Up is displayed in Figure Wlh . It is not a good 
representation of disturbance produced by for example, a flare or eruption, however its 
analysis is useful to understand more realistic excitations. In Figure [T2b the trapped 
energy for the m = 1 mode is plotted. We see that it decreases quite rapidly with tq 
while it smoothly increases with the width of the initial perturbation a. It is possible 
to find the asymptotic behaviour of the energy for a <C rp, 

e-2|Ae|ro 

E = a. (21) 

Due to the strong dependence on rp in general only a small amount of the energy of 
the initial perturbation is trapped in the loop. For example, for a perturbation located 
at a distance similar to the loop length rg = 40/? (L = 50R and a = R), the trapped 
energy is 9.4 X 10~'^_Bo (Eq being the energy of the initial disturbance). This results 
clearly indicates that the loop is able to trap a very small part of the energy of the 
initial disturbance. 

The method to calculate the energy can be easily extended to azimuthally and 
vertically localised perturbations. For example, we can study an azimuthally localised 
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Fig. 12 (a) Example of kink perturbation in the perpendicular plane to the loop axis. For 
this initial disturbance ro = 15-R and a = 2R. In the vertical direction the dependence is 
sin(fcz). White colour corresponds to positive magnetic pressure perturbations while black 
colour represents negative values. The circular white region in the centre represents the loop, 
(b) Trapped energy as a function of ro and a for the m = 1 mode. For this plot L = 50-R and 
pi/pe = 3. The behaviour of the energy with ro and a is given by equation I I21I I in the limit 
ro > R + a. 



perturbation of the form 



The shape of the perturbation and the energy distribution for the different m's are 
plotted in Figure [13] The trapped energy associated with each m is displayed in Fig- 
ure [iSb. It is clear that the energy decreases with m and that the kink mode (m = 1) 
has the largest energy. The energy trapped by the kink mode is between two and three 
orders of magnitude larger than the energy of the first fluting mode (m = 2). Even 
for m = 1 the trapped energy is quite small and, as expected, is lower than that for 
an excitation with a single m since in that case all the initial energy was in one single 
mode. Qualitatively similar results are found when the width of the perturbation in 
the angular direction, 5, is changed. Obviously, the energy of the modes in the initial 
perturbation is different, but again the trapped energy distribution is dominated by 
the kink mode. 

For a longitudinally localised perturbation we assume the following dependence 



where zq is the longitudinal location of the Gaussian and A is the width of the pertur- 
bation in the longitudinal direction. The sinusoidal dependence on z has been included 
in order to satisfy the line-tying condition. This perturbation is represented in Fig- 
ure 114b . and it corresponds to the case zq — L/i, i.e. a perturbation located at a 
quarter of the loop length. Again the kink mode has the largest energy (see Fig. 114b ). 
These results suggests that longitudinal harmonics are in principle more easily excited 
than azimuthal harmonics (which have energy differences of the order 10'^). 




(22) 




(23) 
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Fig. 13 (a) Example of an aaimuthally localised perturbation. The circular white region 
represents the loop, (b) Trapped energy of each eigenmode. Note that the vertical axis is in 
logarithmic scale. The difference in energy between consecutive modes increases with m. For 
example, the trapped energy for the m = 1 is 4 X 10^ larger than the energy of the m = 2. For 
this plot ro = 40R, a = 4R and 5 = 0.5 {L = bOR). 
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x/R 

Fig. 14 (a) Example of longitudinally localised perturbation. The white area represents the 
loop, (b) Trapped energy of each eigenmode (the vertical axis is in logarithmic scale). Note 
that the fundamental mode has the largest energy. For this case ro = 40il, a = 4i? (L = 50il), 
ZQ = L/4 and A = L/10. 



Although there are estimations of the various types of energy released in real flares, 
it is not cleax what amount of energy is associated with MHD waves. Prom the infor- 
mation of the previous model and the information of the coronal loop oscillations, this 
energy can be estimated. For example from the amplitude of the oscillating loops it 
is easy to estimate that they typically have an energy of 10^^ J. If the perturbation 
is located at the loop length, then the energy of the initial disturbance should be 10^ 
times larger, i.e. of the order 10^^ J. This quantity is, in order of magnitude, similar 
to the estimations of the other types of energy release. Thus, the energy of the oscil- 
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Fig. 15 Plot of Vr at the centre of the loop as a function of time for the numerical simulation 
for the initial perturbation given by equation I I17I I. After a short leaky transient, detailed in 
the inner plot, the loop oscillates with the corresponding quasi-mode. The dashed line is a 
fit of the form exp(— t/Tj). In this simulation, rp = 5-R, w = 2R, I = 0.5-R, pi/ pc = 3 and 
R/L = 0.01. 



lating loops can be potentially used as a seismological tool to determine the amount 
of e nergy release d in re al flares (related to fast MHD waves). It is interesting to note 
that lSallai et al.1 1 20081 ) based on a model of forced oscillations (driven by EIT waves) 
have estimated energies of the same order of magnitude. 



4.3 Role of an inhomogeneous layer at the loop boundary 

The effect of a non-uniform layer on the kink oscillatio n has been exten sively studied in 
the literature. The reader is referred to the review bv lGoossenj ( 20081 ) (and references 



therein) about the process of resonant absorption and to its application to explain the 
reported damping of loop oscillations. 

Here, we just want to illustrate by means of a simple time-dependent solution 
the process of resonant absorption that takes place at the loop boundary. For this 
reason, as in the homogeneous cylinder, we set up an initial disturbance given by 
equation (|17|) and we analyse the evolution of the system. Figure [15] shows the radial 
velocity component, Vr, at the loop centre (r = 0) as a function of time for a loop 
with an inhomogeneous layer of I — 0.5-R. We can clearly identify several phases. There 
are initially two large extrema followed by short-period oscillations which are quickly 
attenuated (see inner plot of Fig. ll5|) . After this short leaky transient the loop oscillates 
with a much longer period. Now the amplitude of this mode is attenuated with time 
(compare with Fig. llOb , for the homogeneous loop, i.e. I — 0). This damping is due to 
the conversion of energy from the global kink mode to the localised Alfvenic modes. By 
fitting an exponential function of the form exp(— t/rj) to the envelope of the signal we 
have numerically determined the damping time, rj. There is quite a good agreement 
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Fig. 16 Plot oivr (continuous line) and vg (dot-dashed line) around the resonant position for 
different times for the numerical simulation of Figure 1151 Thin lines correspond to — > oo 
and thick lines to Rm = 10^. Although the behaviour of vg in the inhomogeneous layer is quite 
different for the two values of Rm, the damping times are the same. The vertical dashed lines 
mark the boundaries of the inhomogeneous layer. 



between the time-dependent results and the quasi-mode (which is basically the resistive 
eigenmode in the limit of small dissipation) . The calculation of the quasi-mode provides 
with valuable informati o n abo ut the behaviour of the system, as was also shown by 
iRuderman and Roberts! ( 20021 ) (they solved the time-dependent problem analytically 
using Laplace transforms). 

To show in detail the behaviour of the solutions in the inhomogeneous layer we 
have plotted the time evolution of the radial and the azimuthal velocity components in 
the range < r/R < 2 (Fig. I16|l . For Rm ^ oo we see that small scales are generated 
in time in the region R — 1/2 < r < R + l/2, the amplitude of the azimuthal component 
being much larger than that of the radial component due to the energy conversion. 
The spatial scales in th e inhomogeneous layer decr ease due to the process of phase 
mixing (see for example 'Hevvaer ts and Priestlll983l ) since the local Alfven frequency 
changes with position. Figure [16] also shows the time evolution for the same situation 
but for Rm ~ 10'^ (thick lines). This is an unrealistically small value but it is illustrative 
since it clearly shows the effect of dissipation. Although the damping time is almost 
the same as for the ideal case, now the spatial scales are much larger since the short 
scales have been smoothed out by the effect of diffusion. However, the damping time 
is independent of the value of the resistivity for Rm large. In this limit, resistivity 
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only determines the behaviour in the inhomogeneous layer but it does not affect the 
conversion of energy and in consequence it does not change the damping time. It is 
clear that resistivity plays a secondary role and that the calculations can be done 
even in ideal MHD {Rm — > oo). This property of resonant absorption is particularly 
attractive sine no anomalous Reynolds numbers are required to explain the damping 
of the oscillations. 



4.4 Nonlinear effects on the kink oscillation 

The previous results for the cylindrical loop were derived using the linearised MHD 
equations. Nevertheless, it is also interesting to study the nonlinear regime. One of the 
consequences of a large amplitude (nonlinear) standing kink oscillation is the generation 
of flows along the tube due to the ponderomotive force. This problem was studied 
bv 'Tcrradas an d OfmanI (12004) in the context of coronal loops oscillations (see also 



iHoUwes 1971, ; .Rankin et al.lll99€ ). Due to the nonlinear terms, an initial disturbance 



produces a force imbalance along the loop which results in an up-flow from its legs. 
The accumulation of mass at the top of the loop (for the fundamental standing kink 
mode) can produce strong density enhancements, which secularly grow with time as 



2 



Jc^^l^^j {cotfcos{2k,z), (24) 

where Vq is the amplitude of the initial disturbance. However, if the plasma beta 
is different from zero the gas pressure gradient becomes dominant and inhibits the 
concentration of mass. The density enhancements are then of the form 

where Us — 2kz Cs, being Cs the sound speed. This expression can be potentially used as 
an indirect way to determine the so und speed in the solar coron a if such enhancements 
are reported in the observations (see lTerradas and Ofmanll2004l . for an example). Thus, 
although the zero-/? approximation is in general justified under coronal conditions, the 
role of the gas pressure is important in nonlinear studies of standing kink oscillations 
(and also in propagating kink body modes) since it avoids the eventual destruction of 
the tube due to the ponderomotive force. 

On the other hand, nonlinearity can significantly change the motion of the plasma 
fiuid at certain locations in the tube. It is important to remark that the kink oscillation 
of a cylindrical loop with a sharp or smooth transition involves strong shear motions at 
the tube boundary (see for example Fig. I16|l . Such motions might be liable to develop 
a Kelvin-Helmholtz instability. The most likely place for the KHI to occur is where 
the velo city is largest and t he magnetic field is perpendicular to the velocity (see for 
example iRankin et al.|[l993l ). For the fundamental standing kink mode this is precisely 
the antinode of the velocity, located at half the loop length from the photosphere. In 
order to understand this problem we show the results of so lving the full nonlinear MHD 
equations in three-dimensions (see 'Terrada s et~al]|2008ah . This is the extension of the 
linear problem studied in Sections 14. II and 14.31 Now the initial perturbation is located 
inside the tube. The evolution of the density (for a sharp transition model) and the 
longitudinal component of the vorticity of a representative case in a weak nonlinear 
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Fig. 17 The top half of each panel displays the density (symmetric respect to j/ = 0) and the 
bottom half the vorticity (antisymmetric respect to y = 0) at z = L/2 (half the loop length) 
for different times. In this simulation the following parameters have been used: L = lOR, 
pi/pg=3, vo = O.lwAi, and a full domain of [-6R, -6R] X [-6R,6R] X [0, lOiJ]. The kink 
oscillation develops a Kelvin-Helmholtz instability at the loop boundary. 

regime is shown in Figure [T7] at different time intervals. The initial perturbation pro- 
duces a lateral displacement of the tube in the a;— direction. As in the linear regime 
the tube starts to oscillate around the initial position, but small length scales quickly 
develop at the boundary (see Fig. 117b ). These small scale structures grow with time 
(see Fig. 117b ) and several rolls form at the loop edge. Around x = ±_R and y — the 
density is almost undisturbed because there are no shear motions at this position. At 
later stages of the evolution (Fig. I17t i) the small spatial scales are still localised at the 
boundary. As a result of the instability the shape of the tube at the boundary has been 
considerably changed. 

On the other hand, if instead of a sharp transition a thick inhomogeneous layer 
between the tube and the environment is included then the motions are characterised 
by strong phase-mixed scales (see Fig. (THlfor the linear results). In this situation the 
instability is still present, but it develops more slowly than in the sharp transition 
case. This is due to the fact that the thicker the layer the later the generation of small 
length scales due to phase mixing, and thus the eventual onset of the instability. This 
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process see ms to be related to the devel opment of the KHI for torsional Alfven waves 
described bv lBrowning and PriestI 1 1984h . Interestingly, for thick layers the attenuation 
of the central part of the tube due to resonant absorption (i.e. the damping rate) is 
not significantly altered by the changes at the boundary due to the shear instability. 

Finally, it must be noted that up to now, there is no evidence of such instability from 
the observation of oscillating loops. It is known that magnetic twist, not included in the 
model, might decrease or even suppress the instability since the presence of a magnetic 
field component along the flow stabilises the KHI. Nevertheless, a tube with very large 
azimuthal magnetic field is subject to the instabilities of the linear pinch. Therefore, the 
azimuthal component of the magnetic field of a stable flux tube in the solar corona is 
probably constrained to be in a speciflc range. The absence of observational evidences 
of the instability might be also due to a lack of spatial resolution of the telescopes. 
This issue about this instability is very promising and needs to be investigated in more 
detail in the future from both the theoretical (including twist in the tube models) and 
the observational point of view. 



4.5 Multi-structure 



Now we leave the single tube model to study the collective motion of several tubes 
again in the linear regime. We start with two parallel homogeneous straight cylinders 
of radius R, length L, and separation between centres d. This is th e equivalent to the 
slab configuration analysed in Section [331 In this two-tube problem iLuna et al.l (|2008l ) 
showed that there are four collective fundamental kink-like trapped modes (see also 
Ivan Doorsselaere et al.|[2008l '). The velocity field is more or less uniform in the interior 
of the loops, and so they move basically as a solid body, while the external velocity 
field has a more complex structure. The four velocity field solutions have a well defined 
symmetry with respect to the y-axis. The symmetric mode has the velocity field inside 
the tubes lying in the a:;-direction and is symmetric with respect to the y-axis. This is 
the Sx mode, where S refers to the symmetry of the velocity field and the subscript x 
refers to the direction of the velocity inside the tube. The same nomenclature is used 
for the other modes. The following mode has the velocity inside the cylinders mainly 
in the a;-direction and it is antisymmetric with respect to the i/-axis, we call this mode 
Ax- Similarly, we have the modes Sy and Ay. The pressure field of the Ax and Sy 
modes is symmetric with respect to the y-axis, while that of the Sx and Ay modes is 
antisymmetric. 

In Figure [T^ the time evolution is shown for a pulse in the Vx component initially 
centred on the right loop (see Fig. fTSb). In Figure [T5b. the pulse reaches the left tube 
and passes through it, the system still being in the transient phase. In Figures [TSb and 
I18H the system oscillates in the stationary phase. The oscillatory amplitude in the left 
loop grows with time in the stationary phase, while the amplitude in the right loop 
decreases in the time interval shown in Figures [T5b and I18t l. The left tube begins its 
motion through the interaction with the right loop, i.e. by a transfer of energy from 
the right loop to the left loop. This process is reversed and repeated periodically: once 
the left loop has gained most of the energy retained by the system, so that the right 
loop is almost at rest, the left tube starts giving away its energy to the right cylinder, 
and so on. We find again the beating phenomenon (see Section [3. 3|l . since the initial 
disturbance excites the Sx and Ax modes with the same amplitude. 
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Fig. 18 Time-evolution of the velocity field (arrows) and total pressure field ( contours), for a 
separation between loops d = 6R (in this model a = R). The two circles mark the positions of 
the loops at t = 0. The panels show different evolution times. In (a) the initial condition over 
the velocity field is represented. In panel (b) the tubes are in the transient phase. In panels 
(c) and (d) the system oscillates in the stationary phase with a superposition of the Sx and 
Ax normal modes. 



The extension of the previous model t o multi-stranded config uration but with 
nonuniform layers has been carried out by IXerradas et al.l l|2008 j ). In this work a 
smooth density transitions between the tubes is allowed. In Figure [T9b the two-dimensional 
distribution of the density (the cross section of the composite loop) is plotted for a par- 
ticular configuration. The loop density has an inhomogeneous distribution with quite 
an irregular cross section and boundary. In this so-called "spaghetti model" the dis- 
tance between the strands is quite small and a strong dynamical interaction between 
them takes place. When an initial perturbation is launched in the system the structure 
basically oscillates in the direction of the perturbation. This kink-like motion of the 
bundle of loop results in a strong mode conversion at the inhomogeneous layers. We 
can see the shear motions due to the resonant coupling in Figure 119b (see the phase- 
mixed scales at the boundary). As in the cylindrical model described in Section 33] the 
energy of the kink-like mode is transferred to the layer, this is clear from the energy 
distribution at a particular given time displayed in Figure [20k . In Figure [20b the av- 
erages of the velocity in three difi'erent regions (marked with the labels a, b and c) are 
represented. We see that the signals are almost the same, they have the same dominant 
period and the same damping time. This is due to the global nature of the motion, and 
the attenuation with time is the consequence of the energy conversion due to resonant 
absorption (this result is completely equivalent to that found in the single cylindrical 
inhomogeneous loop model, see Fig. I15|l . These results clearly demonstrate that the 
mode conversion takes place in quite irregular geometries. Regular magnetic surfaces 
are not necessary for this mechanism to work efficiently as was already anticipated by 
iRudermanI (|2003l ) using an elliptical cross section loop model. This indicates that reso- 
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nant absorption is quite a robust damping mechanism and although we have analysed 
a particular system, the behaviour found in the present equilibrium is expected to be 
quite generic of inhomogeneous plasma configurations. The resonant coupling between 
fast and Alfvon modes can hardly be avoided in a real situation since inhomogeneity is 
certainly present in coronal loops, for these reasons, resonant absorption seems to be 
quite a natural damping mechanism. 
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Fig. 19 (a) Cross section of the density of the multi-stranded loop model. The contour 
lines represent curves of constant Alfven frequency, (they also represent contours of constant 
density) . The thick line corresponds to the Alfven frequency that matches the frequency of the 
global mode, (b) Detail of the velocity field at t = 160 ta . 
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Fig. 20 (a) Energy distribution at a given instant (when most of the energy has been trans- 
ferred to the inhomogeneous layers) . The three square boxes labelled with a, b, and c mark the 
regions where the velocity averages have been computed to determine the damping time, (b) 
Time evolution of the spatially averaged Vy component as a function of time. The continuous, 
dot-dashed, and dotted lines correspond to the (normalised) averages in the square regions 
marked with labels a, h, and c in (a), respectively. The profile of the three signals is almost 
the same since the loop is oscillating with a global mode. The damping per period in this 
configuration is r^/P = 2.3. 
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4.6 Curved configurations 

A curved cylindrical model involves a three-dimensional model. Kink oscillations in 
three-dimensional configurations are discussed in this volume by L. Ofman, so the 
reader is referred to this review and references therein. Here we just want to comment 
a few points about these works. One of the striking findings of the three-dimensional 
studies is that, as in some two-dimensional works, the damping rate found in the 
simulations is much stronger that the reported by the observations. This damping is 
usually attributed to the wave leakage due to curvature. However, up to now there 
are no comprehensive numerical studies about this process in 3D. Moreover, if there is 
leakage in three-dimensions, because the frequency of the kink mode is above the local 
Alfven frequency, then it is unavoidable to have an external resonance (basically where 
the frequency of the kink mod e is equal to the local Alfven frequency). This external 
resonance (see lTerradas et all 2006) will also contribute to the attenuation of the kink 
oscillation due to resonant absorption. The implications of this mechanism have not 
been yet addressed in three-dimensional studies. 

5 Conclusions and Discussion 

We have described the main efforts that have been done so far to understand standing 
kink oscillations in coronal loops. We have reviewed the time-dependent studies in 
simple slabs and cylindrical loop models. The effect of curvature and inhomogeneity 
has been also described in some detail, and preliminary studies about the possible 
effects of nonlinearity have been explained. However, the time-dependent modelling of 
kink oscillations is still in a primitive stage and there are several issues that need to 
be addressed, being the attenuation of the oscillations one of the most interesting. 

In particular, the analysis of the damping by wave leakage in two and three- 
dimensional models needs to be investigated in detail, this might help to provide clues 
about the problem of the strong damping found in most of the numerical simulations. 
Nevertheless, here it is worth to draw the attention about the overlooked problem of 
numerical dissipation. Very often, in many numerical studies the effect on the results of 
the intrinsic numerical viscosity of the numerical scheme is not studied or quantified. 
Moreover, in some numerical studies dissipation is explicitly added to stabilise the code. 
The issue of the dissipation is very relevant, specially if we want to extract quantitative 
information about the damping rates with the goal of identifying the possible damping 
mechanism that produce the attenuation. Too much artificial dissipation can simple 
mask the physics behind the damped oscillations and lead to the wrong interpretation 
of the results. 

Another remark about the numerical simulations is that in general in most of the 
works the full set of nonlinear MHD equations are numerically solved. Nevertheless, 
in the majority of the published papers about kink oscillations, the linear regime is 
considered (the amplitude of the waves are very small in comparison with the local 
Alfven or sound speed). Before any problem is attempted to be studied we should 
first ask whether it is worth to solve the full nonlinear equations instead of the simple 
linearised system of equations. One of the advantages of the linear regime is that (for 
simple equilibrium configurations) we can make Fourier analysis in some direction to 
simplify the problem (since that direction is incorporated in the equations through a 
certain wavenumber). This is not possible when the full nonlinear MHD equations are 
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solved. Obviously, if we are interested in nonlinear processes, such as shocks or some 
types of instabilities then it is required to solve the nonlinear MHD equations. 

Finally, we would like to point out that in other branches of sciences, for example 
in computational fluid dynamics, it is a standard practice to attempt to repeat experi- 
ments soon after they are published. This is a necessary aspect of the scientific method 
and it can viewed as an excellent learning tool. We think that it might be i nteresting 
to apply the idea of reproducible research (see the paper about this topic of iLeVeauel 
[2OO6,) in numerical MHD. This would mean to make available the codes used to create 
the figures and tables in a paper in such a way that the reader can download the codes 
and run them to reproduce the results shown in the paper. 

Acknowledgements This paper was inspired by a ISSI workshop held in Bern, Switzerland, 
August 2008. 
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